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We present the results of a time-dependent Dalitz plot analysis of — > tt'^-k~tt'^ decays based 
on a 414fb~^ data sample that contains 449 x 10^ BB pairs. The data were collected on the T(4S') 
resonance with the Belle detector at the KEKB asymmetric energy e'^e~ collider. Combining our 
analysis with information on charged B decay modes, we perform a full Dalitz and isospin analysis 
and obtain a constraint on the quark mixing angle 02, 68° < (j)^ < 95° at the 68.3% confidence level 
for the cf>2 solution consistent with the standard model (SM). A large SM-disfavored region also 
remains. The branching fractions for the decay processes — > p*(770)7r^ and B" p'^(770)7r'^ 
are measured to be (22.6 ± l.l[stat.] ± 4.4[syst.]) x 10"'^ and (3.0 ± 0.5[stat.] ± 0.7[syst.]) x 10"^, 
respectively. These are the first branching fraction measurements of the process B° p(770)7r with 
the lowest resonance p(770) explicitly separated from the radial excitations. 



PACS numbers: 11.30.Er, 12.15.Hh, 13.25.Hw 

1. INTRODUCTION 

In the standard model (SM), CP violation arises from 
an irreducible complex phase in the Cabibbo-Kobayashi- 
Maskawa (CKM) matrix [1]. The SM predicts that mea- 
surement of a time-dependent CP asymmetry between 
the decay rates of B'^ and gives access to the CP vi- 
olating phase in the CKM matrix [2-4]. The angle 02 of 
the CKM unitarity triangle can be measured via the tree 
diagram contribution in b uud decay processes, such 
as B° TT+TT", B° /0^7r=F, or B° p+p' [5]. In these 
decay processes, however, contributions from so-called 
b —I- d penguin diagrams could contaminate the measure- 
ment of (f>2- Snyder and Quinn pointed out that a Dalitz 
plot analysis of B'~' pir, which includes i?° /O^tt^, 
P~T^'^, and B^ — p°7r", offers a unique way to 
determine (f>2 without ambiguity. The Dalitz plot analy- 
sis takes into account a possible contamination from the 
penguin contribution [6]. In addition, an isospin analy- 
sis [7, 8] involving the charged decay modes. i?+ p+Tr" 
and B^ p^"^^ , provides further improvement of the 02 
determination. 

The Belle [9] and BaBar [10] Collaborations recently 
reported the first measurements employing a time- 
dependent Dalitz plot analysis technique. In this paper 
we describe the details of the time-dependent Dalitz plot 
analysis with the Belle detector at the KEKB asymmet- 
ric energy e+e" collider reported in Ref. [9]. We also 
present the first measurements of the branching fractions 
of B" ^ /9±(770)7r=F and B° /9°(770)7r" decay pro- 
cesses obtained from the Dalitz plot analysis, where the 
p{770) is separated from radial excitations. These re- 
sults can be compared with the branching fraction of the 
process B± ^ p"{770)tt^ [11]. 



1-1. KEKB and Belle Detector 

KEKB [12] operates at the T(45) resonance {^/s ~ 
10.58 GeV) with a peak luminosity that exceeds 1.6 x 
10^4 cm-^s-i. At KEKB, the T(45) is produced with 
a Lorentz boost of fS-y = 0.425 nearly along the electron 
beamline (z). Since the B° and B'^ mesons are approxi- 
mately at rest in the T(4S') center-of-mass system (cms). 
At can be determined from the displacement in the z di- 
rection, Az, between the vertices of the two B mesons: 
At ~ Az/Pjc. 

The Belle detector is a large-solid-angle magnetic spec- 
trometer that consists of a silicon vertex detector (SVD), 
a 50-layer central drift chamber (CDC), an array of aero- 
gel threshold Cerenkov counters (ACC), a barrel- like ar- 
rangement of time-of- flight scintillation counters (TOE), 
and an electromagnetic calorimeter comprised of CsI(Tl) 
crystals (ECL) located inside a super-conducting solenoid 
coil that provides a 1.5 T magnetic field. An iron flux- 
return located outside of the coil is instrumented to de- 
tect Kl mesons and to identify muons (KLM). The de- 
tector is described in detail elsewhere [13]. Two inner 
detector configurations were used. A 2.0 cm beampipe 
and a 3-layer silicon vertex detector were used for the 
first data sample of 152 x 10^ BB pairs (DS-I), while a 
1.5 cm beampipe, a 4-layer silicon detector and a small- 
cell inner drift chamber were used to record the remaining 
297 X 10^ BB pairs (DS-II) [14]. 

1-2. Outline of the analysis 

The analysis proceeds in the following steps. First, 
we extract the signal fraction (Sec. 2). We then de- 
termine the sizes and phases of the contributions from 
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radial excitations (Sec. 3). Using the parameters deter- 
mined in the steps above, we perform a time-dependent 
DaUtz plot analysis (Sees. 4-6). The fit results are in- 
terpreted as quasi- two-body CP violation parameters 
(Sec. 7) and as branching fractions of B*^ p^{770)n^ 
and _B" /9*'(770)7r° decays (Sec. 8). We subsequently 
use these results to constrain the CKM angle 02 (Sec. 9). 



1-3. Differential decay width of time-dependent 
Dalitz plot 

We measure the decay process i?" — ^ 7r+7r~7r", where 
we denote the four-momenta of the 7r+, tt", and 7r° by 
p+, and pq, respectively. The invariant-mass squared 
of their combinations 



So = (p+ +P-f 
satisfies the following equation 

s+ + s- + So = "Jro + 2"^^+ + nT'lo 



(1) 



(2) 



by energy and momentum conservation. The differential 
(time-integrated) decay width with respect to the vari- 
ables above (Dalitz plot) is 



dV 



1 



A 



3tt\ 



(27r)3 8m|o 



(3) 



where A 3^ is the Lorentz-invariant amplitude of the 
B°(B0) ^ tt+tt-ttO decay. _ 

In the decay chain T(4S') — > _B"i3" ^ /i.f2, where one 
of the i3's decays into final state /i at time ti and the 
other decays into another final state /2 at time t2, the 
time-dependent amplitude is 



-(r/2+'iA/)(ti+t2 



X \ cos[Amd(ti - t2)/2] {A1A2 - A1A2 



~ism[Amd{h - t2)/2] ( -A1A2 - -A1A2 
\1 P , 



(4) 



Here, p and q define the mass eigenstates of neutral B 
mesons as pB^ ± gS*', with average mass M and width 
r, and mass difference Am^. The width difference is 
assumed to be zero. 

The decay amplitudes are defined as follows, 





= A{B" - 


- fi) , 


(5) 


Al 


= A{W - 
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= A{B° - 
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A2 


= A{W - 


-/2). 


(8) 



In this analysis, we take A^^^ as Ai and choose /2 to be 
a flavor eigenstate, i.e., ^2 = or A2 = 0. Here we call 



the B decaying into /i = fcp ~ 7r+7r~7r° the CP side B 
while the other B is the tag side B, /2 = /tag (/tag)- The 
differential decay width dependence on time difference 
At = tcp — ttag is then. 



dT 



-^l^*l{(|A3.P + |33.|2) 

(|A3^|2- |^3,|2)cos(AmdAt) 
- gtag ■ 2Im (^^A*^A3^^ sin(AmdAt)} dAt 



(9) 



where we assume \q/p\ = 1 {CP and CPT conservation 
in mixing) and \A{B^ ^ /tag)| = \A{W 7,,g)|, and 
integrate over igum = ^cp + ^tag- Here qtag is the fe-fiavor 
charge and qtag = +1 (^1) when the tag-side B decays 
as a i?° (S") fiavor eigenstate. 

Combining the Dalitz plot decay width (3) and the 
time dependent decay width (9), we obtain the time de- 
pendent Dalitz plot decay width 

dr - |^(At;s+,s_)p dAtds+ ds^ , (10) 

where 

|A(Ai;s+,s_)p = e-r|^*l{(|^3.P + P3.|') 

- gtag • (I^SttP - COs(ATOrfAt) (21) 

+ gtag • 2Im (^^Al^A-i^^ sin(AmrfAt)| , 

^37r = ^37r(S-H, S_) , ^3^ = ^3^ (s+ , S_ ) . (12) 



We assume that the tt+tt^tt" decay is domi- 

nated by the i?'' {P'^T amplitudes: P'^Tr~, 
B^ — > p^7r+, and B'^ — > p°7r°, where p can be p(770), 
/9(1450), or p(1700). Although there could exist con- 
tributions from B^ decays into non-pTr 7r"*"7r~7r° final 
states, such as /o(980)7r°, /o(600)7r°, W7r°, and non- 
resonant tt+tt^tt", we confirm that these contributions 
are small; their effects are taken into account as system- 
atic uncertainties (Sec. 6-5). The Dalitz plot amplitude 
A3^(s-i_, s_) can then be written as 

^3.(s+,s-) - Us+,s-)A^ , (13) 

«=(+,-, 0) 

^l3.(s+,s_) = J2 7ni-S+, s^)A\ (14) 

K = ( + ,-,0) 



where A'^(A'^) are complex amplitudes corresponding to 
B"(B°) p+TT- ,p-TT+ ,p°iT° for K = -h,-,0 and the 
functions fK,{s+,S-) incorporate the kinematic and dy- 
namical properties of the B^ decay into a vector p and a 
pseudoscalar tt. The goal of this analysis is to measure 
the complex amplitudes A~^, A~ , A~^, A~ . and A^\ 
we then constrain the CKM angle 02 using these ampli- 
tudes. 
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In contrast to a quasi- two-body CP violation analy- 
sis, the timc-depcndcnt Dalitz analysis includes measure- 
ments of the sizes of the interferences among the final 
states p'^TT", /9^7r+ and p'^tt", and CP- violating asym- 
metries in the mixed final states. In principle, these mea- 
surements allow us to determine all the relative sizes and 
phases of the amplitudes A'^ and A'^, which are related 
to 02 through an isospin relation [7, 8] by 




TT 



A++A^ + 2A" 
A+ + A-+ 2A0 



(15) 



FIG. 1: The relation between three pions in the rest frame of 



Consequently, in the limit of high statistics, we can con- 
strain (p2 without discrete ambiguities. 



1-5. Fitting parameters 



1-4. Kinematics of B° (pTr)° 



The function 
parts as 



can be factorized into two 



(16) 



where FttIsk) and Tj correspond to the lineshape of the 
p and the helicity distribution of the p, respectively. Here 
we assume that a single unique functional form for the 

lineshape F^(s) can be used for all six [15]. Since 
this assumption has no good theoretical or experimental 
foundation, wc check the validity of the assumption with 
data and assign systematic errors. 

The lineshape is parameterized with Brcit-Wigncr 
functions corresponding to the p(770), p(1450), and 
p(1700) resonances: 



-7-i3M^p^^7oo) - (17) 



where the amplitudes /? and 7 (denoting the relative size 
of two resonances) are complex numbers. We use the 
Gounaris-Sakurai (GS) model [16] for the Brcit-Wigncr 
shape of each resonance [17]. 

In the case of a pscudoscalar- vector ( J = 1) decay, Tj 
is given by 



|cos6i> 



(18) 



After (16) is inserted into expressions (13) and (14), 
the coefhcients of Eq. (11) become 



+ 2 E (Rc[/./;] Ut^^ - lm[U;] 
K<CTe{+,-,o} 



+ E (Rc[/./;]/i"; + im[/./;]e) , 

fi<cre{ + ,-,0} 



with 



= 

I. = 

±,Ro (Im) _ 



rRo (Im) 



Imp^A"*] , 

Re (Im) [A" A" ± A" A"*] , 
Re (Im) [A''A''*-{+)A''A''*] 



(20) 



(21) 



(22) 
(23) 
(24) 
(25) 



The 27 coefficients (22)-(25) are the parameters deter- 
mined by the fit [18]. The parameters (22)-(23) and 
(24)-(25) are called noninterfering and interfering pa- 
rameters, respectively. This parameterization allows us 
to describe the differential decay width as a linear com- 
bination of independent functions, whose coefficients are 
the fit parameters in a well-behaved fit. We fix the over- 
all normalization by requiring = 1. Thus, 26 of the 
27 coefficients are free parameters in the fit. 



T+ = -4|p+||p_|cos0+-\ 
Tf = -4|po I |p+ 1 cos , (19) 

ri" = -4|p_||po|cos0-° / 

where pj , p^. arc the three momenta of the and tt'^ in 
the rest frame of p'' (or the tt'tt-' system), and 6^^{= 0^) 
is the angle between pj and pk (see Fig. 1). 



1-6. Square Dalitz plot (SDP) 

The signal and the continuum background e+e^ 
QQ (q = u,d,s,c), which is the dominant background in 
this analysis, populate the kinematic boundaries of the 
usual Dalitz plot as shown in Figs. 2 and 3 (left). Since we 
model part of the Dalitz plot probability density function 
(PDF) with a binned histogram, the part of the distri- 
bution that is concentrated in a narrow region near the 
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FIG. 2; Distribution of signal Monte Carlo (MC) (without de- 
tector efficiency and smearing) in the usual Dalitz plot. The 
dashed line is the kinematic boundary while the hatched re- 
gion corresponds to the region rejected by the mass cut de- 
scribed in Sec. 2. 




FIG. 4: Distribution of signal MC (without detector effi- 
ciency and smearing) in the square Dalitz plot. The solid, 
dashed, and dotted lines correspond to the iso-contours of 
^/sZ = 0.5 GeV, I.OGeV, and 1.5 GeV, respectively, for each 
p'" resonance. 



edge of the usual Dalitz plot is not easy to treat. We 
therefore apply the transformation 



ds^ds^ I det Jldm' d6' , 



(26) 



which defines the square Dalitz plot (SDP) [11]. The new 
coordinates are 



1 

— arccos 2 
TT V m, 



max 




Tinin 



(27) 



given by 
|det J| 

= MP+\\Po\mo ■ 



— — TT sin(7rm') • tt sin(7r0') 



(29) 



where p+ and po are the three momenta of Tr"*" and 7r° in 
the TT+TT" rest frame. Figures 4 and 3 (right) show the 
distributions of the signal and continuum events, respec- 
tively, in the square Dalitz plot. 



1 

TT 



(28) 



where toq = ^/so, rrv^^^ = m — m-^a and to™'" = 2to^+ 
are the kinematic limits of mo, and J is the Jacobian of 
the transformation. The determinant of the Jacobian is 



0.8 ; 
0.6 ; 
0.4 ; 
0.2 ; 
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FIG. 3: Distribution of qq background (from the data Mbc 
sideband) in the usual Dalitz plot (left) and square Dalitz 
plot (right). 



EVENT SELECTION AND 
RECONSTRUCTION 



To reconstruct candidate B" 



'7r° decays. 



charged tracks reconstructed with the CDC and SVD 
are required to originate from the interaction point (IP) 
and to have transverse momenta greater than 0.1 GeV/c. 
Using kaon identification (KID) information, we distin- 
guish charged kaons from pions based on a kaon (pion) 
likehhood Ck(tt) derived from the TOF, ACC and dE/dx 
measurements in the CDC. Tracks that are positively 
identified as electrons are rejected. 

Photons are identified as isolated ECL clusters that 
are not matched to any charged track. We reconstruct 
TT*^ candidates from pairs of photons detected in the barrel 
(end-cap) ECL with E-^ > 0.05 (0.1) GeV, where E^ is 
the photon energy measured with the ECL. Photon pairs 
with momenta greater than 0.1 GeV/c in the laboratory 
frame and with an invariant mass between 0.1178 GeV/c^ 
and 0.1502 GeV/c^, roughly corresponding to ±3(t in the 
mass resolution, are used as tt^ candidates. 
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We identify B meson decays using the energy difference 
AE = E'^^^ — -E^beLm ^'^'^ beam-energy-constraincd 
mass Mbc = Vi^heiJ^ " {pT")^, where E^'^^^ is the 
beam energy in the cms, and E'j^^ and p']^^ are the cms 
energy and momentum, respectively, of the reconstructed 
B candidate. 

We select candidates in a fit region defined as 
-0.2 GeV < A£; < 0.2 GeV and 5.2 GeV/c^ < Mbc < 
5.3GeV/c^. The fit region consists of a signal region 
defined as -0.1 GeV < AE < 0.08 GeV and Mbc > 
5.27GeV/c^, and its complement, called the sideband 
region, which is dominated by background events. 

The vertex position for the B^ 7r+7r~7r° decay 
is reconstructed using charged tracks that have enough 
SVD hits [19]. The /tag vertex is obtained with well- 
reconstructed tracks that are not assigned to fcp- A 
constraint on the interaction-region profile in the plane 
perpendicular to the beam axis is also used with the se- 
lected tracks. 

The 6-flavor of the accompanying B meson is identi- 
fied from inclusive properties of particles that are not 
associated with the reconstructed fcp decay. We 

use two parameters, the 6-flavor charge gtag and r, to 
represent the tagging information [20]. The parameter 
r is an event-by-event, Monte Carlo (MC) determined 
flavor-tagging dilution factor that ranges from r = for 
no flavor discrimination to ?' = 1 for unambiguous fla- 
vor assignment. It is used only to sort data into six r 
intervals. The wrong tag fractions for the six r intervals, 
wi = 1, 2, • • • , 6), and the differences between B^ and 
decays, Awi, are determined using a high-statistics 
control sample of semileptonic and hadronic b c de- 
cays [20-22]. 

The dominant background for the B'^ —f tt^tt^tt'^ sig- 
nal is from continuum. To distinguish these jet-like 
events from the spherical B signal events, we combine 
a set of variables that characterize the event topology 
into a signal (background) likelihood variable £sig (bkg) i 
and impose requirements on the likelihood ratio TZ = 
>Csig/('Csig + -Cbkg)- Due to a correlation between TZ and 
r, these requirements depend on the quality of flavor tag- 
ging. 

When more than one candidate in the same event is 
found in the fit region, we select the best candidate based 
on the reconstructed 7r° mass and TZ. About 30% of the 
signal events have multiple candidates. 

After the best candidate selection, we reconstruct the 
Dalitz variables s+, sq and s_ from 1) the four momenta 
of the 7r+ and 7r~, 2) the helicity angle of the p° (i.e., the 
helicity angle of the tt+tt" system), and 3) the relation of 
Eq. (2). Note that the energy of the 7r° is not explicitly 
used here, which improves the resolution of the Dalitz 
plot variables. We reject candidates that are located in 
one of the following regions in the Dalitz plot: ^/sq > 0.95 
GeV/c^ and ^ > 1.0 GeV/c^ and ^sT > 1.0 GeV/c^; 
^ < 0.55 GeV/c^ or ^ < 0.55 GeV/c^ or ^/sT < 
0.55 GeV/c^ (see Fig. 2). In these regions, the fraction 
of pir signal is small. Moreover, radial excitations 



(the p(1450) and p(1700)) are the dominant contributions 
to B° 7r+7r"7r° in the region with y/s > 1.0 GeV/c^, 
where s is either .s_|-. s_, or sq. Since the amplitudes 
of the radial excitations are in general independent of 
the amplitude of the p{770), they are considered to be 
background in our analysis; vetoing the high mass region 
considerably reduces the systematic uncertainties due to 
their contributions. 

Figure 5 shows the Mbc and AE distributions for 
the reconstructed B^ — > 7r+7r~7r° candidates within the 
AE and Mbc signal regions, respectively. The signal 
yield is determined from an unbinned four-dimensional 
extended-maximum-likelihood fit to the AiJ-Mbc and 
Dalitz plot distribution in the fit region defined above; 
the Dalitz plot distribution is only used for the events 
inside the Ai?-Mbc signal region. The AE-Mhc distribu- 
tion of signal is modeled with binned histograms obtained 
from MC, where the correlation between AE and Mbc, 
the dependence on Pt^o , and the difference between data 
and MC are taken into account. We also take into ac- 
count incorrectly reconstructed signal events, which we 
call self- cross- feed (SCF) and constitutes ~ 20% of the 
signal. In a SCF event, either one of the three pions in 
fcp is swapped with a pion in /tag, or else the n'^ in 
fcp is misreconstructed. We give the details of the AE- 
Mhc and Dalitz plot PDF's of the SCF component in 
appendix A. For continuum, we use the ARGUS param- 
eterization [23] for Mbc and a linear function for AE. The 
AE-Mhc distribution of BB background is modeled by 
binned histograms based on MC. The Dalitz plot distri- 
butions for all components are modeled in the same way 
as the time-dependent fit described later, but integrated 
over the proper time difference, At, and summed over 
the flavor of the tag side B, gtag- The fit yields 971 ± 42 
7r"'"7r~7r'' events in the signal region, where the 
error is statistical only. 




FIG. 5: The Mbc (a) and AE (b) distributions within the 
AE and Mbc signal regions. The histograms are cumulative. 
Solid, dot-dashed, dotted and dashed hatched histograms cor- 
respond to correctly reconstructed signal, SCF, BB, and con- 
tinuum PDFs, respectively. 
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FIG. 6: Mass distributions and fitted lineshapes in p'^Tv~ (left), p~Tv^ (middle), and p'^yr" (right) enhanced regions. The 
histograms are cumulative. Solid, dot-dashed, dotted and dashed hatched histograms correspond to correctly reconstructed 
signal, SCF, BB, and continuum PDFs, respectively. Note that there are feed-downs from other quasi-two-body components 
than those of interest, especially in the high mass regions. For example, the high mass region (mo > l.OGeV/c^) of the p^yr'' 
enhanced region (right) includes large contributions from p^n^ . 



3. DETERMINATION OF THE 
CONTRIBUTIONS FROM RADIAL 
EXCITATIONS 

Although the contributions from radial excitations are 
suppressed by the selections in the Dalitz plot described 
in the previous section, there are still significant contri- 
butions from the long tails of the radial excitations and 
their interferences. We thus need to determine the sizes 
of the radial excitations and their uncertainties to prop- 
erly model the signal PDF's and systematic uncertainties 
associated with their degrees of freedom. 

Using the same data sample as described above but 
performing a time-integrated Dalitz plot fit with a 
wider Dalitz plot acceptance, 0.55GeV/c^ < ^/sq < 
1.5GeV/c2 or ^ < l.SGeV/c^ or ^ < l.SGeV/c^, 
we determine the p lineshape, i.e., the phases and am- 
plitudes of the coefficients /? and 7 in Eq. (17). We use 
these for all of the decay amplitudes. In this fit, we use 
the PDG values [24] for the masses and widths of the 
p(1450) and p(1700). The fit yields 

1/31 = 0.31^0;°^ arg/3=(219l}«)° , 
|7| = 0.08t°:°^ arg7=(l02ti)° . 

The mass distributions and fit results are shown in Fig. 6. 
Figure 7 schematically shows how the radial excitations 
contribute to our fit result. Note that the above values 
are quantities used for time-dependent Dalitz fit and we 
do not regard them as our measurements of /3 and 7. 
This is because these parameters are determined from 
the region where p+7r~ and p~7r+ modes, etc. interfere, 
and they depend on the unfounded common lineshape 
assumption of Eq. (16); hence wc do not give their sys- 
tematic errors. Because statistics are low, we cannot de- 
termine P and 7 for each decay mode without imposing 
the common lineshape assumption. However, we include 
the effect of possible decay-mode dependent differences 
in the values of P and 7 in the systematic errors, which 
are described in Sec. 6-1. 

Thus, it is important to determine the common or av- 
erage lineshape as well as to obtain an upper limit on the 



deviation from the average lineshape for each of the six 
decay amplitudes. For this purpose, we put constraints 
on additional amplitudes that describe 1) the excess in 
the high mass region, y/s > 0.9GeV/c^, where s is ei- 
ther s+, s_, or So; and 2) interferences between radial 
excitations and the lowest resonance p(770) (e.g., inter- 
ferences between p(770)+7r~ and p(1450)~7r+, etc.). The 
nominal fit is performed with the average lineshape de- 
termined above, fixing all of the additional amplitudes 
to zero. When floating the additional amplitudes for the 
other resonances, we obtain results consistent with zero 
for all of the additional amplitudes but with large un- 
certainties compared to the errors for the average line- 
shape parameters above. We use the fit result with the 
additional lineshape parameters floated including their 
uncertainties in the systematic error study. 

4. TIME-DEPENDENT DALITZ PLOT 
ANALYSIS 

To determine the 26 Dalitz plot parameters, we define 
the following event-by-event PDF: 

P(f) = /sigPsig(a') + fsB'PBBi^^) + fm'Pmi^) . (31) 

where Psig, T-'bb ^^"-^ ^99 ^'^'^ PDF's for signal, BB back- 
ground and continuum background, respectively, and 
/sig, f BB ^'^^ ^^"^ corresponding fractions that 
satisfy 

/sig + Ibb + = 1 ■ (32) 

The vector a;, the arguments of the PDF's, corresponds 
to a set of event-by-event variables: 

x= {AE,Mbc;m',e';At,qt^^,l;p^o) . (33) 

A detailed description of the PDF can be found in Ap- 
pendix A. 

With the PDF defined above, we form the likelihood 
function 

C = l[P{x,), (34) 
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where i is an index over events. We maximize C to deter- 
mine the 26 DaUtz plot parameters using the UkeUhood 
function with the signal fraction and the lineshape pa- 
rameters obtained in Sec. 2 and Sec. 3, respectively. 



5. FIT RESULT 

An unbinned maximum likelihood fit to the 2824 events 
in the signal region yields the result listed in Table I. The 
correlation matrix for the 26 parameters, after combining 
statistical and systematic errors, is shown in appendix B. 
Figure 8 shows the projections of the square Dalitz plot 
in data with the fit result superimposed. We also show 
the mass and helicity distribution for each pn enhanced 
region along with projections of the fit (Fig. 9). We find 
that is 4.8 CT above zero, corresponding to clear evi- 
dence for the presence of the decay — > p'^ir^ in agree- 
ment with our previous measurement [25] (see Sec. 8). 
Figure 10 shows the At distributions and background- 
subtracted asymmetries. We define the asymmetry in 
each At bin by {N+ - iV_ )/ {N+ + N_), where (_) cor- 
responds to the background-subtracted number of events 
with (/tag = +1 The p~tt~^ enhanced region shows a 

significant cosine-like asymmetry, arising from a nonzero 
value of UZ . Note that this is not a CP- violating effect, 
since p^n^ is not a CP eigenstate. No sinelike asym- 
metry is observed in any of the p+Tr^, p~7r+ or p^n'^ 
enhanced regions. 

As a check of our fit, we perform the time-dependent 



TABLE I: Results of the time-dependent Dalitz fit. 
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FIG. 7: A schematic figure of the fit result of the lineshape 
and the contributions from radial excitations. Note that our 
definition of -F^(s) does not include the factor l/(l + /3-b7) 
as in Eq. (17). The inset shows the high mass region, 
Wtttt > 1.15GeV/c^, on a semi-log scale where the interfer- 
ence between the p(770) and radial excitations is visible. One 
can see that the p(770) and p(1450) destructively interfere 
with each other near ^/s = m-nw = 1.4 GeV/c^ which means 
that the p(1450) has a large impact on the phase of F-n{s) 
although the absolute value of |J7r(s)| is not much affected. 



Dalitz plot fit with the 5° lifetime floated as a free pa- 
rameter. We obtain 1.41 ± 0.07 ps for the lifetime, where 
the error is statistical only, while the changes of the other 
parameters are very small compared to their statistical 
errors. The lifetime we obtain is consistent with world 
average [24] and thus validates our understanding of the 
PDF's and the background fraction. 




0.2 0.4 0.6 0.8 1 

e' 



0.2 0.4 0.6 0.8 1 

m' 



FIG. 8: Distributions of 6' (a) and m' (b) with fit results. 
The histograms are cumulative. Solid, dot-dashed, dotted 
and dashed hatched histograms correspond to correctly re- 
constructed signal, SCF, BB, and continuum PDFs, respec- 
tively. 
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FIG. 9: Mass (upper) and helicity (lower) distribution of p^-K~ (left), p~Tv^ (middle), and p^n^ (right) enhanced regions. The 
histograms are cumulative. Solid, dot-dashed, dotted and dashed hatched histograms correspond to correctly reconstructed 
signal, SCF, BB, and continuum PDFs, respectively. 




FIG. 10: Proper time distributions of good tag (r > 0.5) regions for /tag = B (upper) and /tag = B (middle upper), in 
p^Tv~ (left), p^tt"*" (middle), p^Tr'^ (right) enhanced regions, where solid (red), dotted, and dashed curves correspond to signal, 
continuum, and BB PDFs. The middle lower and lower plots show the background-subtracted asymmetries in the good tag 
(r > 0.5) and poor tag (r < 0.5) regions, respectively. The significant asymmetry in the p~7r+ enhanced region (middle) 
corresponds to a non-zero value of UZ . 



5-1. Treatment of statistical errors 



With a MC study, we check the pull distributions, 
where the pull is defined as the residual divided by the 
MINOS error. Here, the MINOS error, which corre- 
sponds to the deviation from the best fit parameter when 
— 2 ln(£/£inax) increases by one, is an estimate of the 
statistical error. Although the pull is expected to follow 
a Gaussian distribution with unit width, wc find that 



the width of the pull distribution tends to be signifi- 
cantly larger than one for the interfering parameters due 
to small statistics. We verify with MC that with high 
statistics the widths of the pull distributions are unity. 
To restore the pull width to unity, we multiply the MI- 
NOS errors of the interfering parameters by a factor of 
1.17, which is the average pull width for the interfering 
parameters obtained above, and quote the results as the 
statistical errors. For the non-intcrfcring terms, wc quote 
the MINOS errors without any correction factor. 
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6. SYSTEMATIC UNCERTAINTIES 

Tables II-IV list the systematic errors for the 26 time- 
dependent Dalitz plot parameters. The total systematic 
error is obtained by adding each source of systematic 
uncertainty in quadrature. 



6-1. Radial excitations (p' and p") 

The largest contribution for the interfering parame- 
ters tends to come from radial excitations. The sys- 
tematic error related to the radial excitations (p(1450) 
and p(1700), or p' and p") can be categorized into three 
classes: 1) uncertainties from the lineshape variation, i.e., 
the lineshape difference between each decay amplitude, 2) 
uncertainties in external parameters, ?Tip(]^45o), ^p(u5o)^ 
™p(i700): rp(i7oo), and 3) uncertainties in the common 
lineshape parameters (3 and 7 used for the nominal fit. 

In our nominal fit, we assume all six decay amplitudes 
have the same contribution from (O(1450) and p(1700), 
i.e., we assume Eq. (16). This assumption, however, is 
not well grounded. In general, the contributions from 
p(1450) and p(1700) can be different for each of the decay 
amplitudes and thus the systematic uncertainty from this 
assumption must be addressed. Without the assumption 
about the higher resonances, Eq. (16) becomes 

%^T^'f\{s^), (35) 

where 

+ (/3 + a'^L)W,(i45o)(s) (36) 
+ (7 + A<7^)BWp(i7oo)(s) • 

The variation of the contributions from radial excitations 
is described by non-zero A (3^, and A '7)-, which are 12 
complex variables. We generate various toy MC sam- 
ples, where the input A'^ and A'^ are fixed but the values 
of A and A* 7k are varied randomly according to the 

constraints on A /3„ and A'7k; these constraints arc ob- 
tained from the results in Sec. 3, which are combined with 
the isospin relation [7, 8] to improve the constraints. The 
statistics for each pseudo-experiment are set to be large 
enough so that the statistical uncertainty is negligible. 
We assign the variations and the biases of the fit results 

due to the A (3^, and A'7^ variation as systematic errors. 

For the masses and widths of the /3(1450) and p(1700), 
we use the values from the PDG [24]. To estimate the 
systematic error originating from uncertainties in their 
parameters, we generate toy MC samples varying the 
input masses and widths. We fit them using parame- 
tcrizations with the masses and widths of the nominal 
fit. Here, we vary the masses by twice the PDG error 
(±50MeV/c2 for the p(1450) and ±40MeV/c2 for the 



p(1700)) since the variations between independent exper- 
iments are much larger than the la PDG errors, while we 
vary the widths by the ±1ct PDG errors. We quote the 
mean shift of the toy MC ensemble as the systematic er- 
ror. We also take into account the systematic errors from 
the uncertainties in /3 and 7 for the nominal fit (Eq. (30)) 
in the same way. 

6-2. SCF 

Systematic errors due to SCF are dominated by the un- 
certainty in the difference between data and MC; these 
errors are determined from B — > D^*^ p control samples 
that contain a single 7r° in the final state. We vary the 
amount of SCF by its la error, which is ±100% for the 
CR SCF and tlT°^% for the NR SCF in DS-I (DS- 
II), where DS-I and DS-II denote the subsets of data 
taken with the different detector configurations defined in 
Sec. 1-1. Here, NR represents the neutral-pion-replaced 
SCF and CR represents the charged-pion-replaced SCF. 
We quote the differences from the nominal fit as the sys- 
tematic error. The event fraction for each r region (for 
CR and NR), the wrong tag fractions (for CR) and life- 
time used in the At PDF (for CR), which arc obtained 
from MC, are also varied and the differences in the fit 
results are assigned as a systematic error. 

6-3. Signal Dalitz PDF 

Systematic errors due to the Dalitz PDF for signal 
is mainly from the Dalitz-plot-dependent efficiency. We 
take account of MC statistics in the efficiency and un- 
certainty in the tt'' momentum dependent efficiency cor- 
rection, e'(p^o), obtained from the control samples of 
the decay modes W -> p- D^*'^+ , W -k- D^*^+ , 
B- p-D^*^° and B- tt~D'^*^^. The Dalitz plot efii- 
ciency obtained from MC is found to have a small charge 
asymmetry {'^ 3% at most). We use this asymmetric 
efficiency for our nominal fit. To estimate the system- 
atic error from the asymmetry, we fit the data using a 
symmetric efficiency and conservatively quote twice the 
difference between symmetric and asymmetric efficien- 
cies as the systematic error. The Dalitz plot efficiency is 
r-region dependent and obtained as a product with the 
event fraction in the corresponding region, J~^:^^- €^ irn! ^ ^'), 
using MC. The difference in the fraction for data and MC 
is estimated to be ~ 10% using the B^ D*~tt^ control 
sample. The fractions are varied by ±10% to estimate the 
systematic error. 

6-4. Background Dalitz PDF 

The Dalitz plot for continuum background has an un- 
certainty due to the limited statistics of the sideband 
events, which we use to model the PDF. We estimate the 
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TABLE II: Table of systematic errors (1). The notation "< 0.01" means that the value is small and less than 0.01, and 
thus not visible for the number of significant digits shown here. We calculate the total systematic error including these small 
contributions. 











rr+.He 


rr+.He 








p and p" 


0.01 


0.01 


0.31 


0.19 


0.19 


0.21 


0.39 


0.30 


SCF 


0.01 


0.02 


0.31 


0.09 


0.11 


0.11 


0.12 


0.11 


Signal Dalitz 


0.06 


0.01 


0.15 


0.20 


0.18 


0.13 


0.10 


0.10 


BG Dalitz 


0.02 


0.01 


0.17 


0.11 


0.11 


0.14 


0.10 


0.19 


Other TTTTTT 


0.04 


0.02 


0.06 


0.08 


0.07 


0.10 


0.09 


0.07 


BG fraction 


0.02 


0.01 


0.08 


0.04 


0.07 


0.06 


0.03 


0.11 


Physics 


0.02 


< 0.01 


0.01 


0.01 


0.02 


0.01 


0.01 


0.01 


BG At 


< 0.01 


< 0.01 


0.03 


0.01 


0.01 


0.02 


0.01 


0.01 


Vcrtcxing 


0.03 


0.01 


0.03 


0.05 


0.02 


0.09 


0.05 


0.07 


Resolution 


< 0.01 


< 0.01 


0.03 


0.05 


0.02 


0.02 


0.02 


0.03 


Flavor tagging 


< 0.01 


< 0.01 


< 0.01 


< 0.01 


< 0.01 


< 0.01 


< 0.01 


0.01 


Fit bias 


0.02 


0.02 


0.10 


0.11 


0.07 


0.06 


0.24 


0.22 


TSI 


< 0.01 


< 0.01 


0.01 


0.02 


0.02 


0.02 


0.01 


< 0.01 


Total 


0.09 


0.04 


0.52 


0.35 


0.33 


0.34 


0.51 


0.47 



TABLE III: Table of systematic errors (2). The notation "< 0.01" means that the value is small and less than 0.01, and 
thus not visible for the number of significant digits shown here. We calculate the total systematic error including these small 
contributions. 





^+ 


uz 








I-/-, He 




TT-,lm 




p' and p" 


0.01 


0.02 


0.04 


0.53 


0.29 


0.42 


0.70 


0.31 


0.59 


SCF 


0.02 


0.02 


0.02 


0.09 


0.17 


0.17 


0.13 


0.09 


0.18 


Signal Dalitz 


0.01 


0.02 


0.01 


0.27 


0.20 


0.14 


0.30 


0.15 


0.19 


BG Dalitz 


0.04 


0.03 


0.02 


0.28 


0.32 


0.22 


0.30 


0.20 


0.30 


Other TTTTTT 


0.03 


0.03 


0.02 


0.07 


0.08 


0.12 


0.13 


0.08 


0.08 


BG fraction 


0.02 


0.04 


0.01 


0.18 


0.17 


0.14 


0.22 


0.13 


0.11 


Physics 


0.01 


0.01 


<0.01 


0.03 


0.03 


0.03 


0.04 


0.01 


0.04 


BG At 


<0.01 


<0.01 


<0.01 


0.02 


0.03 


0.02 


0.03 


0.02 


0.03 


Vertexing 


0.02 


0.01 


0.05 


0.18 


0.20 


0.17 


0.08 


0.07 


0.11 


Resolution 


0.01 


0.01 


<0.01 


0.10 


0.14 


0.28 


0.07 


0.11 


0.26 


Flavor tagging 


0.01 


0.01 


<0.01 


0.03 


0.03 


0.03 


0.05 


0.03 


0.02 


Fit bias 


<0.01 


0.02 


<0.01 


0.03 


0.09 


0.02 


0.27 


0.08 


0.26 


TSI 


0.03 


0.03 


0.01 


0.06 


0.03 


0.01 


0.05 


0.04 


0.02 


Total 


0.07 


0.08 


0.08 


0.72 


0.60 


0.65 


0.91 


0.47 


0.83 



TABLE IV: Table of systematic errors (3). The notation "< 0.01" means that the value is small and less than 0.01, and 
thus not visible for the number of significant digits shown here. We calculate the total systematic error including these small 

contributions. 





i+ 


i- 


i(, 






J -1) 






,1m 


p' and p" 


0.02 


0.02 


0.03 


0.82 


0.64 


0.55 


0.46 


0.56 


0.48 


SCF 


0.01 


0.01 


0.01 


0.18 


0.27 


0.10 


0.38 


0.17 


0.14 


Signal Dalitz 


0.01 


0.01 


0.01 


0.28 


0.22 


0.14 


0.27 


0.21 


0.30 


BG Dalitz 


0.01 


0.01 


0.01 


0.29 


0.35 


0.26 


0.28 


0.26 


0.34 


Other TTTTTT 


0.02 


0.03 


0.01 


0.13 


0.10 


0.10 


0.10 


0.13 


0.14 


BG fraction 


0.01 


0.01 


0.01 


0.13 


0.24 


0.19 


0.16 


0.15 


0.25 


Physics 


0.01 


0.01 


<0.01 


0.04 


0.05 


0.03 


0.04 


0.03 


0.05 


BG At 


<0.01 


<0.01 


<0.01 


0.05 


0.04 


0.03 


0.05 


0.04 


0.09 


Vcrtcxing 


0.02 


0.01 


0.03 


0.11 


0.24 


0.09 


0.31 


0.36 


0.16 


Resolution 


0.01 


0.01 


0.01 


0.19 


0.22 


0.15 


0.28 


0.20 


0.23 


Flavor tagging 


<0.01 


<0.01 


<0.01 


0.04 


0.07 


0.04 


0.04 


0.07 


0.03 


Fit bias 


<0.01 


0.01 


<0.01 


0.11 


0.10 


0.41 


0.25 


0.13 


0.18 


TSI 


<0.01 


<0.01 


<0.01 


0.09 


0.04 


0.06 


0.05 


0.18 


0.05 


Total 


0.01 


0.01 


0.05 


0.98 


0.92 


0.80 


0.89 


0.8.5 


0.81 
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uncertainty by performing a toy MC study of sideband 
events. With each MC pseudo-experiment, we model the 
PDF in the same way as wc do for real data. Using the 
PDF, we fit the data in the signal region and quote the 
variation of fit results as the systematic error. The flavor- 
asymmetry parameters for the continuum background, 
which are fitted from sideband events, are varied by their 
uncertainties. Systematic uncertainty from the statistics 
of the BB MC, which is used to model the BB Dalitz 
plot PDF, is also taken into account. 



6-5. B° — > Tr~^iv tt" processes other than 

The primary contribution to the systematic errors of 
the non-interfering parameters tends to come from the 
B'^ 7r+7r^7r° decay processes that are not i?" 
(/97r)'^. We take account of the contributions from 
/o(980)7rO, 50 _^ /o(600)7rO, ^0 ^ ^^0^ ^^^^j 
resonant i?" — > tt'^tt~tt'^. Upper limits on their contribu- 
tions are determined from data, except for — > wtt*^, 
for which we use world averages for B{B'^ ojtt") [26] 
and B{uj tt+tt") [24]. For the mass and width pa- 
rameters of the /o(600) resonance, we use recent mea- 
surements by BES [27], CLEO [28], and E791 [29] and 
take the largest variation. We find no significant signals 
for any of the above decay modes. Using the la up- 
per limits as input, we generate toy MC for each mode 
with the interference between the i?° (P"")" and the 
other B^ — > n^7T~7T° mode taken into account. We ob- 
tain the systematic error by fitting the toy MC assum- 
ing B° [fmf only in the PDF. Within the physi- 
cally allowed regions, we vary the CP violation param- 
eters of the other i?° 7r+7r~7r" modes and the rela- 
tive phase difference between B^ {P'^T ^^-nd the other 
_B" 7r+7r~7r" modes, and use the largest deviation as 
the systematic error for each decay mode. 

In the above procedure, we use relativistic Breit- 
Wigners for the /o(980) and /o(600). To validate the es- 
timated systematic uncertainties, wc investigate possible 
model dependence as follows. For the /o(980), we per- 
form the same procedure using a coupled-channel Breit- 
Wigner [30], which takes account of the opening of the 
KK decay channel, instead of a simple relativistic Breit- 
Wigner. We observe no systematic increase in the un- 
certainties. By changing the model for the /o(980), the 
total systematic error for the "7r+7r~7r° other than pn" 
category increases by at most 10% in the non-interfering 
parameters, which are the only parameters in which the 
systematic error contribution from this category is sig- 
nificant. Thus, we conclude the model dependence of the 
/o(980) resonance parameterization is negligibly small. 
For the /o(600), the situation is more comphcated be- 
cause there are not only possible variations of the reso- 
nance mass spectrum but also uncertainties in the low 
mass 7r~7r+ S'-wave component; this contribution may 
not be modeled by a simple scalar resonance such as the 



/o(600) but by a more sophisticated description known 
as the isT-matrix [31, 32]. To address this issue in a model 
independent and conservative way, we perform the sys- 
tematic error study assuming a hypothetical scalar res- 
onance that has exactly the same mass spectrum as the 
/9"(770). Since this resonance, as a tt+tt" S'-wave contri- 
bution, is maximally similar to the p'^Tr" signal, this pro- 
cedure will lead to systematic uncertainties larger than, 
or at least comparable to, any other model of the tt+tt" 
S'-wave contribution. We find that including the hypo- 
thetical resonance leads to no significant increase of the 
systematic uncertainties compared to those we assigned 
to the /o(600) in our nominal systematic errors. The in- 
crease compared to the case of the /o(600) is at most 30% 
in the non-interfering parameters, which corresponds to 
only a 10% increase of the total systematic error for the 
"7r+7r~7r° other than p7r" category. This is because the 
discrimination power of the p'^ir^ signal from these low 
mass TT+TT" S-wave contribution in general mainly comes 
from the helicity distribution and the result is not sensi- 
tive to the details of the models of their mass distribution. 
Thus, the systematic uncertainties that we have assigned 
for a possible contribution from /o(600)7r° are rea- 

sonable estimates of the systematic uncertainties from 
possible B'^ (low mass tt+tt^ S-wave)7r° contributions 
in general. 

6-6. Background fraction 

Systematic errors due to the event-by-event AE-M\^c 
background fractions arc studied by varying the PDF 
shape parameters; the fraction of continuum background; 
and a correction factor to the signal PDF shape, which 
takes account of the data-MC difference, by ztlcr. We 
also vary the fractions of the BB background, which are 
estimated with MC, by ±50% (±20%) for 6 ^ c (& ^ u) 
processes. 

6-7. Physics parameters 

We use world averages [24, 26] for the following physics 
parameters: Tgo and Am^j (used for signal and BB back- 
ground At), the CKM angles (pi and 02 (used in BB 
background) , and the branching fractions of 5 — > u decay 
modes (used in BB background). The systematic error 
is assigned by varying these parameters by ±1(7. The 
charge asymmetry of B'^ a^n^, for which we use zero 
in the nominal fit, is varied over the physically allowed 
region, i.e., ±1. 

6-8. Background At PDF 

Systematic errors from uncertainties in the background 
At shapes for both continuum and BB backgrounds are 
estimated by varying each parameter by ±lcr. 
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6-9. Vertex reconstruction 

To determine the systematic error that arises from un- 
certainties in the vertex reconstruction, the track and 
vertex selection criteria are varied to search for possi- 
ble systematic biases. In addition to the tracks, the IP 
constraint is also used in the vertex reconstruction with 
the smearing due to B- flight distance taken into account. 
The systematic error due to the IP constraint is estimated 
by varying the smearing by ±10 fim. 



6-10. Resolution function for the At PDF 

Systematic errors due to uncertainties in the resolution 
function are estimated by varying each resolution param- 
eter obtained from data (MC) by ±lcr (±2cr). Systematic 
errors due to uncertainties in the wrong tag fractions are 
also studied by varying the wrong tag fraction individu- 
ally for each r region. 



6-11. Fit bias 

We observed fit biases due to small statistics for some 
of the fitted parameters. Since these biases are much 
smaller than the statistical errors, we do not correct for 
them but rather take them into account in the systematic 
errors. For each parameter, we estimate the size of the 
fit bias from a toy MC study and quote the bias in the 
systematic errors. We also confirm that the bias is con- 
sistent between toy MC and full detector MC simulation. 



6-12. Tag-side interference 

Finally, we investigate the effects of tag-side inter- 
ference (TSI), which is the interference between CKM- 
favored and CKM-supprcssed B ^ D transitions in the 
/tag final state [33]. A small correction to the PDF 
for the signal distribution arises from the interference. 
We estimate the size of the correction using a i?° ^ 
D*~i^v control sample. We then generate MC pseudo- 
experiments and make an ensemble test to obtain the 
systematic biases. 



The time-dependent partial width for the quasi-two- 
body decay process of p^ir^ is given by [34] 



dV 



1 - gtag(C ± AC) cos(AmdAt) (37) 
+ gtag(5 ± A5) sin(AmrfAt)l , 



where the upper (lower) sign is taken for B^ —> p+7r~ 
(p^7r+). The parameters A^J", C, AC, S, and AiS char- 
acterize CP- violating and charge-asymmetric properties 
of 5° — > P^TT^; is a time- and flavor-integrated 

charge asymmetry, C is a flavor-dependent direct CP- 
violation parameter, 5 is a mixing-induced CP-violation 
parameter, and AC and AS are CP-conserving param- 
eters (i.e., non-zero AC or AS does not imply CP vio- 
lation). They are related to the parameters obtained in 
the time-dependent Dalitz plot analysis as 



C 
S 



where 



We obtain 





ir - 

aCP + 


ul 
ul ' 




_ c+ 




_c+ -c- 

2 


(38) 


_ 5+ 




s+ -s- 

2 




c- 






2/_ 
(39) 


A^^ 


= -0.12 ±0 


05 ±0.04, 


(40) 


c 


= -0.13 ±0 


09 ±0.05 , 


(41) 


AC 


= +0.36 ±0 


10 ±0.05 , 


(42) 


s 


= +0.06 ±0 


13 ±0.05 , 


(43) 


AS 


= -0.08 ±0 


13 ±0.05 , 


(44) 



where first and second errors are statistical and system- 
atic, respectively. The correlation matrix is shown in 
Table V. 

The angle 02^, which equals (/)2 in the no-penguin limit, 
can be defined as [35] 



.off 
*2 ^ 



off'.- 



(45) 



7. QUASI- TWO-BODY PARAMETERS 

One can deduce quasi-two-body CP- violation parame- 
ters from the fit result of the time-dependent Dalitz plot 
analysis. In this section, we obtain the quasi-two-body 
CP- violation parameters from the results for the U and 
/ parameters determined in the previous section. 



TABLE V: Correlation matrix of the quasi-two-body param- 
eters, with statistical and systematic errors combined. 







C 


AC 


5 


A5 




+1.00 










C 


-0.17 


+1.00 








AC 


+0.09 


+0.16 


+1.00 






5 


+0.01 


-0.02 


-0.00 


+1.00 




AS 


-0.00 


-0.01 


-0.02 


+0.29 


+1.00 
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with 



2<^f ± 6 



S±AS 



arcsm 



v/l-(C±AC)2 



and 



5 = arg {A-*A+) 
Our measurement gives 



(88.0 ±3.9 ±1.7)' 



(46) 



(47) 



(48) 



where (j)2^ 



90° would give values of S and AS con- 
sistent with zero. There also exists a mirror solution 
= (2.0 ± 3.9 ± 1.7)° due to the two-fold ambiguity in 
the arcsine. In addition, other solutions (/)|^ ~ 45° and 
135° are also allowed in principle. The additional solu- 
tions correspond to cases where 2(f>2^''^ + 5 and 2(/)2*'~ — S 
differ by ~ 180°; they can be excluded by including a 
weak theoretical assumption (flavor SU(3) or QCD fac- 
torization implies a much smaller value) [35]. The mea- 
sured (f>2^ can be used to constrain 02 in a model de- 
pendent way, using a theoretical assumption that puts a 
limit on the difference <j)2 — 4>2^ [35, 36]. 

The direct CP-violation parameters for the process 
p'^TT^, Ap~ and -4p^, are defined as 



p-7r+) -r(B" -> p+n^) 



and 



p-TT+)+T{BO p+TT-) 

p+n-) - T{B" -> p-'K+) 



p+TT-) + r(BO p-7r^ 



(49) 



(50) 



One can transform the parameters A'^^ , C, and AC into 
the direct CP violation parameters as 



A-1 



A 



CP 



C 



1 + AC 



(51) 
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FIG. 11: Contour plot of the confidence level for the direct 
CP violation paxameters .4p^^ vs. ApJ'. 



where Ap0.^o and SpO^o are the parameters to be mea- 
sured. They are calculated from the parameters fitted in 
the time-dependent Dalitz plot analysis as 



We obtain 



ApO^^O — - 



ApO-^0 = 
SpO^O = 



2Io 



-0.49 ±0.36 ±0.28 
-0.17 ±0.57 ±0.35 



(56) 



(57) 
(58) 



with a correlation coefficient of —0.08. We observe a 
very small correlation between the quasi-two-body CP- 
violation parameters of the processes P" p^n^ and 
P" p°7r", whose absolute values are less than about 
0.02. Our measurement of ApO^o is consistent with the 
previous measuremeirt from Belle [25]. 



A-+ = 



^,^/-C-^^/AC 



AC 



We obtain 



A^ 
A 



pTT 



-0.21 ±0.08 ±0.04 
-0.08 ±0.16 ±0.11 



(52) 



(53) 
(54) 



with a correlation coefficient of +0.47. Our result differs 
from the case of no direct CP asymmetry (yl^^ = and 



A- 



0) by 2.3 standard deviations (Fig. 11). 



We also measure the CP violating parameters of the 
quasi-two-body P° p°7r° decay process. The time- 
dependent partial width for the process is given as 



dV 
dAt 



|At|/r^ 



1 + gtag.4p0^o cos(AmdAt) 



QtagSpO^o sin(AmdAt) 



(55) 



8. BRANCHING FRACTION MEASUREMENTS 



The number of P" 
in Sec. 2 is 



TT+TT tt" events, Nsia, obtained 



iVsig = 971 ± 42 



The branching fraction for P° {pT^T 
eluding radial excitations is given as 



i3(P" ^ pTT^") 



]\J _ gDct gVoto g/KID ' 



(59) 



(60) 



where e^®' and 6^®*° are the average detection efficiency 
and the efficiency corresponding to the Dalitz veto (i.e., 
the upper and lower bounds on s+, s_, and sq), respec- 
tively; e'KiD jg -(-jj^j, efficiency correction factor to take 
account of the KID difference between data and MC; 



15 



and iV^-g is the number of BB pairs produced, where 
we assume equal numbers of B^B'^ and B^B~ pairs. 
The branching fraction for the decay to a ground state 
p'^(770)7r°' is given as 



[K,al = [±,T],[0,0] 



( I' 



(61) 



where /p^^ and /p«(77o) are the fractions of i?° Paii"''^ 
among all i?° pTr^" — > tt^tt~tt'^ and that of _B" 
P'^{770)tt'^ among i?° — > p^u^'^j respectively. Here, 
p^jj symbolically represents the total contribution from 
p'=(770), p«(1450), and p'"(1700). The fraction /p~^ 

is calculated from the parameters t/+ and [7^^^°*'^™'' 
(k, (T = +,—,0) in Table I, while the coefhcients rep- 
resenting contributions from radial excitations 7) in 
Eq. (30) determine /p'<(77o)- More details of the formal- 
ism can be found in Appendix C. The values for these 
coefficients are shown in Table VI. Note that /p«(77o) is 
close to 1 and 1 1 — fp^ (779) | is much smaller than the error 
in our result. This means that it is reasonable to com- 
pare the central values of our result with the preceding 
branching fraction measurements that do not separate 
the contribution from radial excitations. However, the 
errors are not directly comparable. 

From (59) and the coefficients in Table VI, we obtain 



B{B° 



pn' 



all\ 



(25.8 ±1.2 ±3.6) X 10" 



(62) 



and 

BiB° p±(770)7rT) = (22.6 ± 1.1 ± 4.4) x 10"^ , 
B{B^ /9°(770)7r") = (3.0 ± 0.5 ± 0.7) x 10"^ . ^^^^ 

Here, the first and second errors correspond to the sta- 
tistical and systematic errors, respectively, where statis- 
tical errors include the contributions from the statisti- 
cal errors on the number of events (Ngig) and uncertain- 
ties in the Dalitz parameters (U^ and J/^^^^^™-*). The 
correlation coefficient for the statistical errors between 
BiB° ^ p±(770)7r=F) and B{B^ ^ /9"(770)7r0) is -0.09. 
The branching fractions obtained are consistent with our 
previous measurements [25, 37]. 



TABLE VI: Summary of the coefficients used in the branching 
fraction measurement. 



^Dct 


0.10 


^Vcto 


0.84 


g/KID 


0.96 


Pall 


0.89 




0.12 


/p± (770) 


0.99 


/p0{770) 


0.99 



8-1. Systematic uncertainties 

Table VII summarizes the systematic errors for the 
branching fraction measurement. We discuss each item 
in the table in the following. 



1. Radial excitations (p' and p" ) 

As in the time-dependent analysis, the systematic un- 
certainty related to the radial excitations come from 1) 
uncertainties due to lineshape variation, 2) uncertainties 
in masses and widths of the p(770) and its radial excita- 
tions, and 3) the common lineshape parameters. Their 
impacts are estimated in the same manner as done for 
the time-dependent analysis. The impact from possible 
lineshape variation, which is constrained by our data, 
is estimated by an MC study. We adopt the same un- 
certainties for the masses and widths of the p{770) and 
radial excitations as in the time-dependent analysis and 
estimate their impact on the branching fractions. The er- 
rors due to common lineshape parameters (/?, 7) are also 
taken into account. 

The uncertainty for this category is sizable since the 
p(770) and radial excitations are separated in the branch- 
ing fraction analysis. This is an essential difference from 
the time-dependent analysis, where we do not separate 
the p(770) and the radial excitations. 



2. Physics parameters 

This category includes the systematic error from the 
uncertainties in the branching fractions of the BB back- 
ground components as well as those due to the possi- 
ble contribution from the B^ tt+tt^tt'' processes other 
than B" ^ (p7r)0. 



3. Fit 

This category includes the uncertainties related to the 
extended unbinned maximum likelihood fit, except for 



TABLE VII: Summary table of the systematic errors for the 
branching fraction measurements. A common factor of x 10~^ 
is omitted for simplicity. 





all 

pn 


p=^(770)7rT 


p"(770)7r" 


p' and p" 


±2.8 


±3.9 


±0.5 


Physics Parameters 


±0.1 


±0.1 


±0.0 


Fit 


±1.7 


±1.5 


±0.2 


Detection Efficiency 


±1.4 


±1.3 


±0.2 


TDPA Systematic 


±0.2 


±0.4 


±0.4 


Number of BB 


±0.3 


±0.3 


±0.0 


Total 


±3.6 


±4.4 


±0.7 
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TABLE VIII: Summary of the systematic uncertainties re- 
lated to the detection efficiency. 

7r° Detection 4.7% 

Kaon Identification (KID) 0.4% 

Continuum Suppression 2.3% 

Vertex Reconstruction 1.8% 

Dalitz Efiiciency 0.8% 



the items included in above two categories. It consists of 
the uncertainties from the modeling of PDF used in the 
fit, where the possible uncertainty in the SCF component 
has a sizable impact; the fraction of the BB background 
component, which is fixed in the nominal fit; and the fit 
bias. 



4- Detection efficiency 

The largest components of the detection efficiency sys- 
tematic uncertainty are the differences between data and 
MC. We consider differences in 7r° detection efficiency 
with t:^ momentum dependence, the KID efficiency cor- 
rection (e'^^^), the continuum event suppression cut, and 
the vertex reconstruction efficiency. All are estimated 
from control sample studies. Small uncertainty due to 
the limited statistics of MC used in calculation of the ef- 
ficiency is also taken into account. Table VIII shows the 
breakdown of the systematic uncertainty contributions 
from the above items. 



5. TDPA systematic errors 

The systematic errors in the Dalitz plot parame- 
ters obtained in the time-dependent Dalitz plot analy- 
sis (TDPA) listed in Table I propagate to the branching 
fractions. 



including (j)2'. 

9 =(6 complex amplitudes = 12d.o.f.) + 02 

— (1 global phase) — (1 global normalization) (64) 

— (1 isospin relation = 2d.o.f.) , 

where we make use of an isospin relation that relates 
neutral B decay processes only [7, 8]. Parameterizing 
the 6 complex amplitudes with 9 free parameters, we 
form a function using the 26 measurements from our 
time-dependent Dalitz plot analysis as constraints. We 
first optimize all the 9 parameters to obtain a minimum 
X^j Xmin! then scan 02 from 0° to 180° optimiz- 
ing the other 8 parameters, whose resultant minima are 
defined as x^(02); the difference Ax^(02) is defined as 
Ax^(02) = X^{.4>2) - Xmin- Performing a toy MC study 
following the procedure described in Ref. [38], we ob- 
tain the 1 — C.L. plot in Fig. 12 (dotted line) from the 
Ax^(02)[39]. 

In addition to the 26 observables obtained from our 
time-dependent Dalitz plot analysis, we use the branch- 
ing fraction B{B'^ pTr^'^) obtained in Sec. 8 and the 
following world average branching fractions and asymme- 
tries: B{B+ p+TfO), A{B+ p+ttO), B{B+ ^ p°7r+), 
and A{B^ p'^t^^) [26], which arc not correlated with 
our 26 observables. With the 31 measurements above, we 
perform a full combined Dalitz and isospin (pentagon) 
analysis. Having 5 related decay modes, we have 12 free 
parameters including 02: 

12 =(10 complex amplitudes — 20d.o.f.) + 02 

— (1 global phase) (65) 

— (4 isospin relations = 8d.o.f.) . 

The detail of the construction can be found in ap- 
pendix D. The Xmin obtained is 10.2, which is reason- 
able for 31 (measurements) — 12 (free parameters) = 19 
degrees of freedom. Following the same procedure as 
above, we obtain the 1 — C.L. plot in Fig. 12 (solid line). 
We obtain 68° < 02 < 95° as the 68.3% confidence in- 
terval consistent with the SM expectation. Several SM- 



6. Number of BB pairs 

The number of accumulated BB pairs and its uncer- 
tainty are N^-^ = (449.3 ±5.7) x 10^, assuming an equal 
production rate for charged and neutral BB pairs from 
the T(4S'). The uncertainty in the number of BB pairs 
propagates to the branching fraction and is taken into 
account as a systematic error. 



9. CONSTRAINT ON THE CKM ANGLE 02 

We constrain the CKM angle 02 from our analysis fol- 
lowing the procedure described in Ref. [6]. With three 
— * decay modes, we have 9 free parameters 




30 60 90 120 150 180 

412 (degrees) 

FIG. 12: 1 — C.L. vs. 02. Dotted and solid curves correspond 
to the result from the time-dependent Dalitz plot analysis only 
and that from the Dalitz and an isospin (pentagon) combined 
analysis, respectively. 
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disfavored region (0° < (j)2 < 5°, 25° < 02 < 32°, and 
108° <(j)2 < 180°) are also allowed. 



10. CONCLUSION 



Using 414 fb ^ of data we have performed a time- 



events (i = NR,CR): 



dependent Dalitz plot analysis of the 



-TT° de- 



cay mode. Combining our analysis and information from 
charged B decay modes, a full Dalitz plot and isospin 
analysis is performed to obtain a constraint on (j)2 in 
a model-independent way. We obtain 68° < 02 < 95° 
at the 68.3% confidence interval for the solution consis- 
tent with the SM expectation. However, a large CKM- 
disfavored region also remains. In principle, with more 
data we may be able to remove all the additional 02 
solutions. From the result of the Dalitz plot analysis, 
we also obtain the branching fractions for the decays 
B° p±(770)7r=F and /9"(770)7r". These are 

the first branching fraction measurements of these pro- 
cesses with the lowest resonance p{77Q) explicitly sepa- 
rated from the radial excitations. 
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APPENDIX A: PDFS FOR TIME-DEPENDENT 
DALITZ PLOT ANALYSIS 

In this section, we describe the details of the PDF for 
each component, which appear in Eq. (31). 

1. Signal PDF 



i=NR,CR 



V^{X) 



true + X]i=NR,CR 



(Al) 



where NR and CR represent the 7r° (neutral) replaced 
and TT^ (charged) replaced SCF's, respectively, and N 
are the integrals of the PDF's. 



a. PDF for correctly reconstructed events 

In terms of the event fractions for the l-ih flavor tag- 
ging region {Tl^^^), the Dalitz plot dependent efficiency 
(e'), the TT*' momentum-dependent efficiency correction 
taking account of the difference between data and MC 
(e'), wrong-tag fractions (lu;), and the differences in 
wrong-tag fractions between B^ and B^ (Au;/), the PDF 
for correctly reconstructed events is given by 

Ptruc(^) =.^truc ' Ptrue ( A£;, Mbc ; P^o ) 

• etrue(m', 9'; I) e'{p^o) ■ |det J(m', 9')\ (A2) 

• Ptiue{m' ,0';At,qte,g;l) , 

where the Dalitz-At PDF 7'truc('7i', 6*'; Ai, gtag; corre- 
sponds to the right-hand side of Eq. (11) and is 

Vuuc{rn',9';At,qtai.;l) 
e-l'^*l/^BO 



(l-<ZtagAlZ;i)(|^3.r + |^3.r) 



- gtag(l - 2wi) (lAa^p - lAg^p) cos(AmdAt) 

+ 9tag(l - 2u;,)2Im (^^Aa^A^^^ sin(AmrfAi)| . 

(A3) 

For the At PDF, the above equation is convolved with 
the resolution function [19]. The terms I^stP ± I^SttP 
and Im (^^A^^^A^^^^ are expanded as in Eqs. (20) and 
(21). The A£'-Afbc PDF is normahzed such that 



dAE dM^, VtrueiAE, M^,-p^o) = 1 (Vp^o) , 
signal region 

.(A4) 

since we define the Dalitz plot efficiency for events inside 
the signal region. Since the PDF's in the At-gtag direc- 
tion are also normalized to be unity, the integral inside 
the signal region, Ntmc, is, 



(A5) 



where 



The PDF for signal events consists of a PDF for the 
correctly reconstructed events T'truo and PDFs for SCF 
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^ dAt dAE dMbc / / dm' d9' Vt_ 

„ J J ^signal region J -/SDP, Veto 



T\,^^ f [ dm' dO' etruc(™', e'- 1) e'ip^o) |dctJ| {{As^f + 

JJSBP, Veto 



(A6) 



r 



and the correlation between p^o and m' is properly taken 
into account in the integration on the last line. The 
notation J /g^p y^^.^ dm' dO' means integration over the 
square Dalitz plot with the vetoed region in the Dalitz 
plot taken into account. 

The tt" momentum dependent Ai?-Mbc PDF, 
Vtruc{AE , Mbc] Pno) , is modeled using MC-simulated 
events in a binned histogram interpolated in the p^o 
direction, to which a small correction obtained with 
— > is applied to account for the difference 

between MC and data. 

The Dalitz plot distribution is smeared and distorted 
by detection efffciencies and detector resolutions. We ob- 
tain the signal Dalitz plot efficiency from MC to take 
the former into account. We introduce a dependence 
of the efficiency on the r region, e[j.^g, since a signifi- 
canct dependence is observed in MC. Small corrections, 
e'(p^o), are also applied to the MC-determined efficiency 
to account for differences between MC and data. We 
use W p-D'-*^+, W -> n-D*-, B- p^D" and 
B~ TT~ decays to obtain the correction factors. 
The smearing in the Dalitz plot due to the finite detec- 
tor resolutions is small compared to the widths of p(770) 
resonances; the smearing is confirmed by MC to be a 
negligibly small effect. 



h. PDF for SCF events 

Approximately 20% of signal candidates are SCFs, 
which are subdivided into - 4% NR SCF and ~ 16% CR 



SCF. It is therefore important to develop a model that 
describes the SCF component well. The time-dependent 
PDF for SCF events is defined as 

r^{AE,Ah,:m',9';At, gtag;0 

=Tl ■ V,{AE, il/bc; s,) ■ r,{m',0'; At, gtag) , (A7) 
{i = NR, CR) 

with 



Viim' ,e'\At,qtag) 
e-|At|M 



4t,; 



[{l~q,,^AwC)Tl''%m',e') 



- gtag(l - 2wi') V^'°'{m', 6') cos(AmdAt) 
+ gtag(l - 2we) Vf'^m', 0') sin(AmrfAO} ■ 

. . 

where J-^ is the event fraction in each tagging r-bin. 

The A£'-A/bc PDF is normalized inside the signal re- 
gion as 



dAEdMy,,V^{AE,M^,■s,) = l (Vs.). (A9) 

signal region 



As the PDF's in the Ai-gtag direction are normalized to 
unity and = 1, the integral inside the signal region, 

A/",, is 



J 



dAt dAEdMbc 

J J signal region 



/ / drn' dB' V\{AE, M^,; m', 6'; At, gtag) = // 

"'JSDP,Voto J Js 



dm' d0' rj-'^^m'^l 

SDP, Veto 



(AlO) 



r 



We find that the AE-Mhc distribution for SCF has a 
sizable correlation with Dalitz plot variables, but only in 
one of its two dimensions. We thus introduce a model 
with dependences on the Dalitz plot variable Si. The 
variable scR = s± = max(s-|_,s_) is used, because the 
CR SCF can be divided into a 7r+ replaced SCF and a 
7r~ replaced SCF, where s„ (s+) is used for tt^ (tt") 
replaced SCF. Here, we exploit the fact that almost all 



of the 7r+ (tt ) replaced SCF distributes in the region of 
s+ < s_ (s+ > s_). For the NR SCF, snr = sq. This 
parameterization models the correlation quite well, with 
each of the parameters Si related to the kinematics of 
replaced tracks. 

Since track (tt) replacement changes the measured 
kinematic variables, the SCF events "migrate" in the 
Dalitz plot from the correct (or generated) position to 
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the observed position. Using MC, we determine resolu- 
tion functions ^i(TOobs' ^obs! "^gen) ^gen) describe this 
"migration" effect, where {m'^bs^9'ohs) and (m'^en^^'gcn) 
are the observed and the generated (correct) positions 
in the DaUtz plot, respectively. The resolution function 
satisfies the normalization condition of 



/ / '^™obs'^^obs -^i("^obs' ^obs' "^gcn'^gcn) ^ ^ 
J JSDP 



(All) 



(ym' 6' ) 



Together with the efficiency function ei{m'g^^,6'g^^), 
which is also obtained with MC, the Dalitz plot PDF 
for SCF is described as 

riim',e') = \iR,-e,)®Pl^ (m',0') 

SDP 

• ^iifTT-gcm 9 gen) ' ^phys('^gcni ^gcn) i 

i =CR, NR , j = Life, Cos, Sin ) 

(A12) 



where 



pLifc / ' f)' ) 

physV gen' gen/ 

pCos / / nl \ 

physV geni gen/ 



|detJ|(|A3,|2 + |^3,|2) 

|detJ|(|A3.p-|33.n 



^phys("^;en'Cn) = |detJ|2Im ( ^A^^A*,, 



(A13) 



For the NR SCF, the shape of the At PDF defined in 
Eq. (A8) is exactly the same as correctly reconstructed 
signal, i.e., tnr = tbo, wi^^ = wi, and Aw;^^ = Aw/, 
since the replaced 7r° is not used for either vertexing or 
flavor tagging. On the other hand, for the CR SCF, the 
At PDF is different from correctly reconstructed signal, 
since the replaced tt^ is used for both vertexing and flavor 
tagging. Thus, we use MC-simulated CR SCF events 
to obtain tcr, wp^, and Aw;^^, which are different 
from those of correctly reconstructed signal events. In 
particular, Aui;*"^^ is opposite in sign for the tt"*" and tt^ 
replaced SCFs, due to the fact that the replaced tends 
to be directly used for flavor tagging in the slow pion 
category. 



Continuum PDF 



The PDF for the continuum background is 



E 

-l + qt^gA\m',9') 



T\-q ■ V\-JyAE, A/be) • n^(m', 0'; AS, Afbc) 



(A14) 



ng(At) , 



where J-\^-^ is the event fraction for each r region obtained 
in the signal yield fit. All the terms on the right hand 
side of the equation are normalized to be unity so that 



I q^. ^ signal region 



signal region 

X ndm' d9' Vqg{AE, Mbc; m', 9'; At, gtag; /) = 1- 

/SDP, Veto 

(A15) 

Since the allowed kinematic region is dependent on AE 
and Afbc, the Dalitz plot distribution is dependent on 
AE and A/bc- We define a A£:-A'fbc independent PDF, 
^9g('^scaic' ^'): whcrc m'^caic a re-dcfincd SDP variable 
with the kinematic effect taken into account as 



scale 



where 



1 

■ arccos 2- 



Too — TO™'" 



-, mm 
'■0 



AE + AA/bc 



AA/bc = A/bc - TO^o 



(A16) 



(A17) 



Using the Ai?-A/bc independent PDF, 
Vqqim', 9'; AE, M\,c) is described as 



Vqq{m',9';AE,Ah,} = 



1 



Nqq{AE + AA/b. 

sin(7rTO') 



sin(7rTO^^^iJ 



1 V ■'Pqq{m'sc^\c,9') 



(A18) 



for the region, m™'" < toq < min(TO™' , ;/(,q 
AAi^bc) {Pqq = otherwise), where A/'5g(A£' + AA^bc) and 
sin(7rTO')/ sin(7rTOgj.g^]g) are a normalization factor and the 
Jacobian for the parameter transformation Tn'^^^-^^ m! , 
respectively. We obtain the Pqqi'tn'^caic^^') distribution 
from data in part of the sideband region. — O.lGeV < 
AE < 0.2 GeV and 5.2GeV/c2_< A/bc < 5.26GeV/c2, 
where the contribution from BB background is negligi- 
ble. 

Since we find significant flavor asymmetry depending 
on the location in the Dalitz plot, we introduce the fol- 
lowing term to take account of it: 



l + gtagA'(m',0') 



(A19) 



which is r region dependent. The asymmetry is anti- 
symmetric in the direction of 9', i.e., A^{m!,9') > 
{A\m',9') < 0) in the region of 9' > 0.5 {9 < 0.5), 
and the size of the asymmetry is at most ~ 20% in the 
best r region. This effect is due to the jet-like topology 
of continuum events; when an event has a high momen- 
tum TT^ (tt^) on the CP side, the highest momentum tt 
on the tag side tends to have -I- (— ) charge. The highest 
momentum tt on the tag side with + (— ) charge tags the 
flavor as i?° (B'^). Since an event with a high momen- 
tum TT~ (7r+) resides in the region 9' > 0.5 {9' < 0.5), a 
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continuum event in the region 6' > 0.5 {9' < 0.5) tends 
to be tagged as We again parameterize the 

asymmetry A'-{m',9') in a Ai?-Mbc independent way as 

A'{m\ 0') = A\m\ 0'; AE, M^,) = A'(m^,,i„ 0') , 

(A20) 

and model it with a two-dimensional polynomial, whose 
coefBcicnts are determined by a fit to data in the AE- 
Mbc sideband region. 



APPENDIX B: CORRELATION MATRIX OF 
THE FIT RESULT 

Tables IX-XI show the correlation matrix for the 26 
parameters determined in the time-dependent Dalitz plot 
analysis, corresponding to the total error matrix with 
statistical and systematic error matrices combined. We 
assume no correlations for the systematic errors. 



3. BB background PDF 

The treatment of BB background is different for CP 
eigenstate modes and flavor-specific or charged modes. 
The PDF for the CP eigenstate modes is 

V'g-^{AE,Ah,;m',e';At,qt,g) 

= ^BB ■ 'PsBi^E, Afbc) • VBBim'j') ■ VBs{At, gtag; 0, 

(A21) 

where Vg-g{At,qta,g) is a time-dependent CP violation 
PDF normalized as 

J dAt VssiAt, gtag; = 1- (VO (A22) 



For the flavor specific or charged modes, the PDF is 
VegiAE, il/bc; m', 0'; At, gtag; 



^BB-'^BBi^E,M^,) 



(A23) 



where the Dalitz plot PDF Vgg{m' ,6';qrcc) is depen- 
dent on the true flavor of the CP (fully reconstructed) 
side, grec, and the time dependent part is a mixing PDF 
(lifetime PDF with flavor asymmetry) for flavor spcciflc 
(charged) modes. The At PDF is normalized as 

5^ ^ y"dAt Vg-siAt, 9tag, grcc; = 1- (VO (A24) 

9tag 91-00 

The AE-M^c PDF and Dalitz plot PDF are obtained 
mode- by- mode from MC. The Dalitz plot PDF of the 
CP eigenstate modes is assumed to have the following 
symmetry 



VB-^{m',6')^Vs-^{m',l-0') 



(A25) 



while that of flavor specific and charged modes is assumed 
to have the following symmetry 

VBB{m, 0': Qr-ec) = VBBim', 1 - 0' ; -qrec) ■ (A26) 

The total PDF of the BB background is a linear com- 
bination of each mode with efficiencies and branching 
fractions taken into account. 



APPENDIX C: FORMALISM FOR BRANCHING 
FRACTION MEASUREMENT 

In this section, we describe the formalism for the 
factors used in the branching fraction calculation of 
Eqs. (60) and (61). 



1. Detection Efficiency 

Since the detection efficiency is Dalitz plot dependent, 
we use the detection efficiency averaged over Dalitz plot, 
e^®', for the branching fraction measurement. From 
Eqs. (Al), (A2), (A3), (A7), (A8), and (A12), the PDF 
integrated over A£'-Mbc (in the signal region) and At, 
and summed over gtag and I is 



Vsieim',0' 



nruc(m',0') + E,=CR, NR^^K.^') 



A/true + Si=CR, NR 



(CI) 



with 

Vtruc{m',0') 

= ;^.Ft',,,etruc(m',0';/)e'(p,o)Pphys(m',0') , 

(C2) 

V^{m', 0') = Y.Tl [(e, • i?,) ® Pphys] (m', 0') , (C3) 
and 



-Pphys (m- ,0 ) 

= |J(m',0')l(l^3.|' + |A3.n 

|2 



J{m',0')\\S2u+\f, 

+ (C/+/<=Re[/./:] - t/+/-Im[/„/: 



(C4) 



Consequently, Vsigim' ,0') can be symbolically rewritten 
as 



Vsisim\0') 



[e"«'^fphys](m^gO 

A/truG + Z]i=CR, NR-^« 



(C5) 



21 



TABLE IX: Correlation matrix (1) of the 26 fitted parameters, with statistical and systematic errors combined. 
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where 2. Dalitz Veto EfRciency 

The efficiency corresponding to the DaUtz veto, gV°*°, 



«> ^'phys] (to', 0') is simply calculated as 

J2Kuc^truc{ni', 9'- Oe'(p.o)Pphys(TO', 9') 



(C6) // dmWPphysK,^') 

X — ^ X — ^ J r in,. Veto ^SDP. Veto / r^r\\ 

^ Y,^l[{^^■R^)®P^''''']{r^^\e'). e =-77 — — • (C9) 

z=CR, NR i dmde PphysiTO , 9 ) 

J JSDP, Whole 



Note that Ri is normalized as Eq. (All). 

Consequently, the detection efficiency averaged over 
the Dalitz plot e^"^* is 



3. Fraction of p^uTt" 



The fraction of Paii^'^ normalized to 



,Det 



dm'd9' 

SDP, Veto 



• P'P^^%m' ,9') 



[[ dm'd9' PP^y%m\9') 

J J SDP, Veto 
A/true + X) j=CR, NR 



dm'd9' PP^y''{m',9') 

SDP, Veto 



(C7) 
(C8) 



dm'd9'P!^{m',9') 



SDP. Whole 



dm'd9' Pphys{m\ 9') 

SDP, Whole 



(CIO) 



where 



P^^:{m',9')^\J{m',9')\U:\f.\^ 



(CU) 
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TABLE X; Correlation matrix (2) of the 26 fitted parameters, with statistical and systematic errors combined. 
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TABLE XL Correlation matrix (3) of the 26 fitted parameters, with statistical and systematic errors combined. 
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4. Fraction of p''(770)7r'^ 

The fraction of p'^{770)t:'' normalized to B° 



fp'^ {770) 



SDP, Whole 



dm'dO' Pl^^^{m',e') 

SDP, Whole 



where 



The function /k'"™'' is defined as 



5. Summary 



By using the expressions described above, Eqs. (60) 
and (61) arc rewritten as 



(C12) 



dm'd0'Pphys(m',0') 

^Vsig J J SDP, Whole 

^BB -'^truc + X]i=CR, NR-^i 



(C15) 



(C13) 



(CM) 



and 



corresponding to the p(770) part of defined in 
Eqs. (16) and (17). 



Bip'' {770)71") 



dm'de' P^'^^.{m',0') 

SDP, Whole 



-^BB -'^true + X]i=CR, NR 



(C16) 



23 



APPENDIX D: METHOD OF 02 CONSTRAINT 
1. Formalism 

We define amplitudes as 



and 





= A{B^ - 


> p+TT ) , 


(Dl) 




= A{B° - 






A" 


= A{B^ - 




(D3) 




= A{B+ - 




(D4) 




= A{B+ - 




(D5) 


A+ 


= ^A{W- 
P 




(D6) 


Ar 


= ^A{W- 
P 




(D7) 


A° 


= ^A{W- 
P 




(D8) 


A-o 


^ UiB- 
P 




(D9) 




^ U{B- 
P 




(DIO) 



These amplitudes are obtained from 1) 26 measurements 
determined in the time-dependent Dalitz plot analysis as 
well as 2) branching fractions and asymmetry measure- 
ments, and give a constraint on (/)2. 

Equations (22)-(25) define the relations between the 
amplitudes for the neutral modes and the parameters 
determined in the time-dependent Dalitz plot analysis. 
The relations between the branching fractions and asym- 
metries, and the amplitudes are 

B{pTT^") = c- {\A^\' + m')-rBO , (Dll) 



The amplitudes are expected to follow SU{2) isospin 
symmetry to a good approximation [7, 8] 



A+ + A- + 2A" = A+ + A- + 2A" 

= V2(A+o + = V2{A-° + I"-) 



(D16) 



A+" - A°+ - V2{A+ -A-)^A-^- A°- - V2{A- ~A+) , 

(D17) 

where 



and I"- = 6-2**^0- 



(D18) 



Note that there is an inconsistency in Eq. (D17) between 
Rcf. [7] and Ref . [8] ; we follow the treatment of Ref. [7] . 



2. Parameterization 

Here we give two examples of the parameterization of 
the amplitudes. The first example may be more intu- 
itive, while the second example is well behaved in the 
fit. The results are independent of the parameterizations 
with respect to the constraint on (f>2- 

In the following, we set the constant c to be unity and 
discard the normalization condition jA+p -I- = 1 

instead. This is equivalent to letting c be a free parameter 
in the (j)2 constraint fit, keeping the normalization 
l^+P = 1. We adopt the former for simplicity. Note that 
Eqs. (22)-(25) become 



I. 

7-r±,Re(Im) 
jRe(Im) 



(|A«|2± /N , 

Im [A'' A"*] /N , 

Re(Im) [A'^A" ±A''A"] /N , 

Re(Im) \A''A''*-{+)A''A''*] /N 



B{p+n')^c-{\A+r + \A-r)-rB^ 



(D12) 



= I 4 + 12 



N=\A 



B{p\+)=c-{\A"+\' + \A"-\^)-rB. , 



A{p+7:'>) 



|^_0|2_ |^+0|2 



1^ 



-012 



|A+"|2 



_ \AO-\^-\AO^? 



|AO-|2 + |ylO+|2 



(D13) 



(D14) 



(D15) 



where c is a constant and the lifetimes Tgo and are 
introduced to take account of the total width difference 
between B° and B~^ . Note that we do not use quasi- 
two-body parameters related to neutral modes except for 
B{p^TT^), since they are included in the Dalitz plot pa- 
rameters. 



in this case. 



a. Amplitude parameterization 



We can parameterize the amplitudes as follows [7] 



A+ 
A- 

A" 



P- 
1 



(D19) 
(D20) 



e-"^^T° - -(P+ + P-) , (D21) 



(D22) 



2T" - T+°) -P+ + P- , 
(D23) 
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and 



V2A-" 



^ e+''t"^T- +P- , 



(D24) 
(D25) 

(D26) 

(D27) 



2r" - 



(D28) 

where the overall phase is fixed with the convention 
ImT+ = 0. Thus, there are 6 complex amplitudes, 
r+,T-,r°,P+,p-, and T+°, corresponding to 11 de- 
grees of freedom; and (/)2, corresponding to 12 degrees of 
freedom in total. This parameterization automatically 
satisfies the isospin relations without loss of generality, 
i.e., the isospin relations are the only assumption here. 



h. Geometric parameterization 




We can parameterize the amplitudes using the geo- 
metric arrangement of Fig. 13 that satisfies the isospin 
relation of Eq. (D16). This figure is equivalent to Fig. 3 
of Ref . [8] , except that the sides corresponding to — > 
/9~7r+ and p^n^ are swapped. This difference is not 

physically significant. Wc apply this modification only to 
obtain a better behaved parameterization; the parame- 
terization here uses the angles a;_ and 0^ related to the 
process B^{B'^) p~i^^ i which are better behaved than 
those related to P"(P°) p°7r°. 

To parameterize the amplitudes, we use 02 and the 
following 11 geometric parameters: 

ti>+, cj-, a;', 0^,9-, 6+, b-, b', a+, a-, L, (D29) 

where b and a imply branching fraction and asymmetry, 
respectively. In terms of these parameters, the ampli- 
tudes can be described as follows 



A+ 


= e'('^+ 




a+)/2, 


(D30) 


A+ 






a+)/2 , 


(D31) 


A~ 




+»-/2)v/5-(l- 


«-)/2, 


(D32) 


A- 


= e''^"" 


-«-/2)v/6-(l4 


-a_)/2. 


(D33) 




= (L- 


A+-A-)/2, 




(D34) 


1° 


= (L- 


I+-I-)/2, 




(D35) 


A°+ 


= e^^' \ 


/&V2, 




(D36) 
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- A"+ , 




(D37) 
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(D38) 



and 



L 

2y/2 



A+" - A 



0+ 



V2(A- 



V2{A+ - 
-I+)l/2 



A- 



(D39) 



FIG. 13: Complex pentagons formed from the _B — > pvr decay 
amplitudes. 



Equation (D39) exploits the isospin relation of Eq. (D17), 
which Fig. 13 does not incorporate geometrically. The 
phase 02 enters when the A's are converted into A's with 
Eq. (D18). When we perform the analysis only with the 
time-dependent Dalitz plot observables and without the 
information from charged decay modes, we remove the 
parameters lo' and 6' from the fit and fix L to be a con- 
stant. 

This geometric parameterization has a substantial ad- 
vantage in terms of required computational resources, 
compared to the parameterization based on the T and 
P amplitudes described in the previous section. In the 
procedure to constrain 02, the minimum has to be 
calculated for each value of 02- To avoid local minima, 
initial values of the parameters in the minimization have 
to be scanned. This inflates the computing time, which 
increases exponentially with the number of parameters. 
However, the number of parameters to be scanned de- 
creases in the geometric parameterization. Among the 
11 parameters except for 02, five of them, 6_, 6', a.)-, 
and a_ , are related to the branching fractions and asym- 
metries. Since in most cases they do not have multiple 
solutions, we do not have to scan their initial values. In 
addition, the optimum initial value for L can also be 
determined using other parameters and 6o, the nominal 
branching fraction of P" p^i^^ , from the following re- 
lation 



6o= L-e''^+v/W2-e'"-v/M2 



(D40) 



up to a two-fold ambiguity. Here &o is calculated using 



25 



the input parameters as 



ur + ui c • tbo 



(D41) 



based on Eqs. (22) and (Dll). The exphcit sohition for 
the optimal initial value of L is 



L ~ Re7 ± ybo — (Im-y) . 
( where 7 = e*'^+ ^/b+/2 + e"^- ^/b-/2 ) 



(D42) 



When 60 — (Im7) < 0, there is no real- valued solution 
and L = Re7 is the optimum initial value. With the op- 
timum values calculated above, the initial value of L does 
not have to be scanned, except for the two-fold ambiguity. 
Consequently, the number of parameters to be scanned 
in this parameterization is only five, corresponding to 
W-)-, W-, w', and while 10 of the 11 parameters 
have to be scanned in the T and P amplitude parame- 
terization. This leads to a substantial reduction of the 
computational resources required. 
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